Equation of state for partially ionized carbon at high temperatures 
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Equation of state for partially ionized carbon at temperatures T > 10 5 K is calculated in a 
wide range of densities, using the method of free energy minimization in the framework of the 
chemical picture of plasmas. The free energy model includes the internal partition functions of 
bound species. The latter are calculated by a self-consistent treatment of each ionization stage in the 
plasma environment taking into account pressure ionization. The long-range Coulomb interactions 
between ions and screening of the ions by free electrons are included using our previously published 
analytical model. 
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I. INTRODUCTION 

The understanding of the physical properties of mat- 
ter at high densities and temperatures is important for 
the fundamental physics as well as for various physical 
and astrophysical applications. Since the 1980s the the- 
oretical interest in matter under such unusual conditions 
(e.g., Refs. P, 0, 0, H) has been enhanced by labora- 
tory developments like high-power short duration lasers, 
shock-induced plasmas, inertial confinement implosions, 
or exploding metal wires (e.g., Refs. 0, IS 01 HI @ ) • In the 
astrophysical domain the calculation of the equation of 
state (EOS) for stellar partial ionization zones is a par- 
ticularly challenging problem. In these zones the elec- 
trons and different ionic species cannot be regarded as 
simple ideal gases: Coulomb interactions, bound-state 
level shifts, pressure ionization, and electron degeneracy 
should be taken into account. In this paper, we calculate 
the EOS for carbon at temperatures 10 5 K < T < 10 7 K 
in a wide range of densities p. Such an EOS is required, 
e.g., for modeling inner envelopes of carbon-rich white 
dwarfs or outer envelopes of neutron stars. 

An EOS calculation in the partial ionization regime is 
not possible without approximations. For astrophysical 
simulations, these approximations should not violate the 
thermodynamic consistency. The free energy minimiza- 
tion method [l2L IT^ | allows one to include the complex 
physics in the model and ensures the consistency. This 
method has the great advantage to identify the various 
contributions to the free energy illustrative of various 
physical effects (see, e.g., Ref. [lj], for a discussion). 

Free-energy models which carefully include the non- 
ideal effects have been proposed for fluid hydrogen [l5l ] 
and helium [l6l Il7|: the EOS tables for these elements, 
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which cover a pressure and temperature range appropri- 
ate for low-mass stars, brown dwarfs, and giant planets 
have been calculated in 

111 

For heavier elements, a 
similarly detailed EOS is lacking. Up to now, the best 
available thermodynamically consistent EOS for carbon 
covering the stellar pressure ionization zones was the 
EOS developed by Fontaine, Graboske, and Van Horn 
in the 1970s 0] (FGV) and amended in 1990s 0. 
We shall call these two versions FGV77 and FGV99, 
respectively. This EOS has been calculated by differ- 
ent methods in different p — T domains. At relatively 
low densities (e.g., p < (0.01 — 1) g cm~ 3 for 10 5 K 
< T < 10 6 K), the ionization equilibrium has been ob- 
tained by the free-energy minimization technique. At 
densities above several grams per cubic centimeter, the 
Thomas-Fermi model has been employed. At interme- 
diate densities, in particular in the various regimes of 
pressure ionization, the EOS was interpolated between 
these two regions. Clearly, the accuracy of the EOS 
in the interpolation region can be called into question. 
Moreover, the Thomas-Fermi model may be inaccurate 
at p < 10 3 g cm -3 , where the pressure is not sufficiently 
high to force the complete ionization of carbon, as we 
shall see below. 

Extension of the free energy minimization technique to 
9 ^ 0-1 § cm" 3 is complicated because of the growing im- 
portance of nonideal contributions to the free energy and 
the onset of pressure ionization. The latter is difficult 
to treat in the framework of the "chemical picture" of 
plasmas, which assumes that different ion species can be 
clearly identified (see, e.g., Refs. 0,H3,0|, for discus- 
sion). On the other hand, EOS calculations within the 
more rigorous "physical picture," quite successful at rela- 
tively low p (e.g., |22(), become prohibitively complicated 
at such high densities. First principle approaches based 
on path integral Monte Carlo (PIMC) or molecular 
dynamics (MD) calculations are computationnaly highly 
expensive. These methods also suffer from some diffi- 
culties. Indeed, the sign or node problem for the PIMC 
method or the use of effective pair potentials for MD 
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simulations restrict their applicability (see however 0|). 
In any case, a comparison with our results will be in- 
structive, but, to the best of our knowledge, no PIMC or 
MD data for carbon in the temperature-density range of 
interest in this paper has been published yet. 

In this paper we present an EOS model which relies 
on the free energy minimization in the framework of the 
chemical picture and extends to arbitrarily high densi- 
ties across the pressure ionization region without inter- 
polation. This allows us to obtain not only the thermo- 
dynamic functions, but also number fractions for every 
ionization stage. We treat the long-range interactions 
in the system of charged particles (ions and electrons) 
using the th eory previously developed for fully ionized 
plasmas [2^, l26j |. The contribution of the internal elec- 
tronic structure of the ions embedded in the dense plasma 
is calculated using a scheme 27] which self-consistently: 
(i) builds separate models for different ionization stages 
in the plasma, taking into account the real structure of 
bound states (configurations, LS terms); (ii) uses Boltz- 
mann statistics to sum up the internal partition functions 
of these ions; (iii) takes into account spreading of bound 
states into energy bands as they are pressure ionized; 
and (iv) treats quantum mechanically the free electron 
background around each ion thus resolving resonances. 
Points (i) and (ii) make our model different from average 
atom ones. The closest "ion-in-plasma" theoretical model 
is that of Pain and Blenski |24j, where ions are treated 
separately (using superconfigurations) , but screening is 
introduced through a Thomas-Fermi approach for the 
free electrons. The applicability of our model is tested 
by numerical calculations of thermodynamic functions, 
which we compare with the FGV models. 

In Sec. [H] we present the free energy model. The tech- 
nique for the calculation of thermodynamic functions at 
equilibrium is described in Sec. lIIII In Sec. II VI we discuss 
the results of the EOS calculations for carbon plasma, 
and in Sec. we give conclusions. 



II. FREE ENERGY MODEL 

Consider a plasma consisting of N e free electrons and 
Ni heavy ions with numbers of bound electrons v from 
to Z nuc (where Z nuc is the element charge number) 
in a volume V. Let us write the total Hclmholtz free 
energy as F tot = F e + F { + F cx , where F he denote 
the ideal free energy of ions and free electrons, respec- 
tively, and F ex is the excess (nonideal) part, which 
arises from interactions. Fi is the free energy of an 
ideal Boltzmann gas mixture, which can be written 
as Fi = Nik B T [ln( ni \f)-l] - S miyi T : where A 4 = 
{2-kK 2 /mik^T) 1 ' 2 is the thermal wavelength of the ions, 
mi is the ion mass, S m \x = —Nik^^2 v x v \nx v is the 
entropy of mixing, and x v — N v /Ni is the number frac- 
tion of the ions of the ^-th type [Y^ u x " ~ ^or * ne 
electrons at arbitrary degeneracy, F e can be expressed 
through Fermi-Dirac integrals and approximated by an- 



alytical formulae |25j . The main complication is the cal- 
culation of the nonideal term, which is quite nontrivial 
at high densities. It includes a contribution due to the 
building of localized bound states of the ions, and the 
long range Coulomb interactions between these ions and 
free electrons. We write 

^cx F ee -\- Fa -\- Fi e -\- F[ n i : (1) 

where the first three terms represent the contributions of 
electron-electron, ion-ion, and ion-electron interactions, 
respectively, and Fi nt is the contribution due to the in- 
ternal degrees of freedom of the ions, that involves sums 
over bound states. Equation does not imply a strict 
separation of the terms on its right-hand side: No strict 
definition of free and bound electrons nor ions exists in a 
dense plasma. In general, the terms must be interdepen- 
dent and evolve in a correlated way. Our approach to this 
difficulty consists in calculating self-consistent models for 
the ions embedded in the plasma, coupling them with a 
model for the long range interaction, and minimizing the 
resulting total free energy F tot . 

A. Free energy of a fully ionized plasma 

A fully ionized electron-ion plasma which contains only 
one ion species is characterized by three parameters: the 
ion charge Ze, the electron density parameter r s , and the 
ion Coulomb coupling parameter T: 

ao' V 3 J k B T 

where n e is the electron number density and ao = 
h 2 /m e e 2 is the Bohr radius. The Hclmholtz free energy of 
the fully ionized plasma is described by analytical fitting 
formulae [25l .26), which are applicable at high densities 
( r s 1j arbitrary T) or high temperatures (small T). 
When neither r s nor T are small, the plasma cannot be 
considered as fully ionized. 

In a multicomponent fully ionized, dense plasma with 
different ion charges Z u e, the "linear mixing rule" has 
been shown to be very accurate [2^, |2^, |3(j : 

F^(Ni,V,T,{x v },{Z v }) = NikaT^XuU, (3) 

where 

/„ = / e fi x (n e , T, Z v ) = i£| Xi/=1 /JVi*BT (4) 

is obtained from = F ee +Fa+Fi e (the superscript "fi" 
indicates full ionization). In Eq. n e takes the value 
implied by the electroneutrality: n e = niZ, where rij = 
Ni/V is the ion number density, and Ze = '^2, v x v Z v e is 
the mean ion charge. An effective ion Coulomb parame- 
ter for a multicomponent plasma is obtained by replacing 

z 5 / 3 with Y.v x ^ lz in E q- @ for r. 
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B. Bound-state contribution to the free energy 

In order to evaluate -Fi n t, we calculate the ionic struc- 
ture in the plasma using the scheme described in ■ It 
is based on the ion-sphere approximation, which replaces 
the actual plasma environment for every ion by the sta- 
tistically averaged boundary conditions for the electron 
wave functions within a spherical volume centered at the 
ionic nucleus. At present we do not include neutral atoms 
(y = Z nuc ), which is justified at the temperatures and 
densities where the ionization degree of the plasma is 
high. For each ion containing v bound electrons, a ra- 
dius of the ion sphere R v is determined self-consistently 
from the requirement that the sphere is overall electri- 
cally neutral. The Hamiltonian for the ion v is written 
as H v = Yh=i M r i) + W v , where 



hu(r) = - Y V 2 + V: t (r) + V?(r), 



(5) 



r, (6) 



i<3 



Vf" is the potential due to the plasma on the ion v, that 
must be determined self-consistently, W v is responsible 
for the LS splitting of spectroscopic terms, and V£ t is a 
scaled Thomas-Fermi potential of the nucleus and v — 1 
bound electrons jU- Note that V£ t disappears in H v . It 
is used to build an effective one-electron Hamiltonian h Vl 
which generates a one-electron wave functions basis. The 
coordinate parts 4>nim °f these functions are obtained 
from the Schrodinger equation 



hv^nlm — e vnl4'nlm( r )- 



(7) 



Then H v is diagonalized in a subspace of Slater deter- 
minants generated by a set of ipnim- The ^-electron 
energies of the bound states are well approximated as 



4nl), Where E L aild 

e tni are calculated for the isolated ion, and a = 
(nl)i(nl)2 ■ ■ ■ {nl) u 2S+1 L defines a particular LS term 
of a configuration. The separation of H v into parts (J5J 
and ijfjj allows one to capture the plasma effects in one- 
electron energies and wave functions through Eq. JSJ, 
while the ^-electron structure is retained through the 
contribution W v . The boundary condition at R v for Eq. 
does not noticeably affect E va except near the den- 
sities where the corresponding term a becomes pressure- 
ionized. The latter case will be addressed below. 

The free electron density rif (r) and the potential (r) 
are determined self-consistently, using the local density 
approximation of the density functional theory. The one- 
electron wave functions tp"i m (r) of the partial scattering 
waves are calculated from the Schrodinger equation 



TO 



+Vr(r) + Kc (r) Vei m {r) = e ^ m (r). (8) 




23 24 25 

log n° (cm -3 ) 

FIG. 1: Neutrality volumes of carbon ions, v v normalized to 
the fiducial number density of free electrons n° (see text), 
as functions of n° for T = 2.32 x 10 5 K (solid lines) and 
T = 3.71 x 10 6 K (dotted lines) for carbon ions with v bound 
electrons {y = 0, 1, ... , 5). The curves for the three lowest 
degrees of ionization (y = 3, 4, 5) end at the n° e values at which 
there remain no relevant bound states (e_ become positive for 
2s state in ls 2 2s and ls 2 2s 2 , and 2p state in ls 2 2s 2 2p). 



Here, V£. is the exchange-correlation potential [33, 
and V f v are obtained from the Poisson equation: V 2 V^ = 
—Airn^e 2 , V 2 V[ = — Airrif e 2 , and the number densities 
n£(r) and rif(r) are calculated as the squared moduli 
of the wave functions for the bound and free electrons, 
respectively, summed with the statistical weights appro- 
priate for a given T. For the bound electrons, these 
weights are proportional to w va d va ey^p{—E va /k^,T), 
where d va = (25 + 1)(2L + 1) is the level degener- 
acy, and w va is an occupation probability defined be- 
low. The density of states per unit volume for the 
Zth partial wave of the free electrons at a given energy 
gvi(e) is determined with account of the contribution 
from resonances (Friedel terms; see Ref. [27| for details). 
The energy distribution of the free electrons is assumed 
oc g v i{e)/(l + cxp[(e — /j, e ) / k^T]) . The free parameters 
of the model are T and the electron chemical potential 
/i e . In thermodynamic equilibrium, [i e is the same for 
all ions, but as different ionization stages have different 
neutrality sphere radii as well as different numbers of 
neutralizing free electrons inside them, fi e can be related 
to the mean free electron density only after the global 
free-energy minimization fSec. lIIljl has given the relative 
populations of the ions. 

The neutrality of the ion sphere is ensured by the self- 
consistent determination of R,. such that 



_d_ 

dr 



+ V h v (r) + Vr(r) 



(9) 



r=R.„ 



Associated with this radius is the neutrality volume 
v v = AttRI/3. In the model of a uniform electron back- 
ground, that neglects the interactions of free electrons 



4 



with ions, one has v v = = (Z nuc — v)/n®, where 
n°(fj, e ,T) = V~ 1 dF e /d/j, e \v,T is the number density of 
free electrons in the uniform gas model. With allowance 
for interactions of the free electrons with the ions and 
bound electrons, v v deviates from v®, as illustrated in 
Fig. H The drops of the plotted curves at certain densi- 
ties, which are especially sharp at the lower temperature, 
are the consequence of pressure ionization of separate lev- 
els: when a nl level of ion v crosses the continuum limit 
and appears as a resonance in the neighboring ionization 
state v — 1, the latter ion sphere shrinks to compensate 
this increase in the free electron density of states. 

With increasing /j, e (or n°), the radius R v decreases, 
the wave functions VnZm( r ) become distorted, and the 
energies e un i spread into a band. We estimate a band 
width by solving Eq. Q with two alternative boundary 
conditions: either ipnim(^) = °' or ®' l l ) nlm/® r = at 
r = R v . These two conditions give two energies which 
we interpret as the upper (e + ) and lower (e_) edges of 
the band (Fig. |2J). Eventually e+ becomes positive. We 
interpret the electrons with < e < e + as quasifree 
and exclude them from the internal partition function 
of the ion. We introduce an occupation probability w un i, 
equal to the statistical weight of electrons with e < (the 
significance and thermodynamic meaning of occupation 
probabilities in the chemi cal p icture of plasmas has been 
discussed, e.g., in. [2(J, |33, |34( ) . Assuming for the bands 
the Hubbard [3!| density of states 



9(e) = ^V5 2 -(e-c) 2 , S = 
ird z 

we obtain, for e_ < < e+, 



(10) 



l y — 

Wvni = V 1 

I 7T 



y 2 arcsin(y), (H a ) 



2/=(e+ + e_)/(e + -e_). 



(lib) 



The occupation probability of a term a is w va = 
ri(n2)6a w vni- For all electron shells nl, except the K 
shell, e_ becomes positive at sufficiently high /i e ; in this 
case w vn i = and the bound state disappears. The low- 
est curves in Fig. ^ end at the densities where the bound 
states cease to exist in the plasma for a given v. 

For the K shell of H- and He-like ions, e_ is negative at 
any density. Asymptotically, at large fi e (small R u ) and 
for a uniform density n e , e_ ~ — 0.3 (Z nuc — v)e 2 /R v , 
while e+ ~ 0.5m~ 1 (irh/R„) 2 . Therefore, at high densi- 
ties w vnl ~ (4\/6/7r 4 )[(2 , lmc - is)R u /5a ] 3/2 oc n^ 1/2 . 

The electrons, that populate the bands overlapping 
with the continuum, are delocalized, and thus add to 
the free-electron degeneracy through the Pauli exclusion 
principle. It means that at a given n.j the presence of 
such electrons increases /i e . Conversely, at a fixed /z e , the 
electrons that are pushed into continuum have a larger 
neutrality volume, corresponding to the unscreened shell. 
Since the share of these quasifree electrons is (1— w V a)i an 
effective ion charge for such a partially delocalized state 
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FIG. 2: Monoelectronic energies of He-like carbon (y = 2) 
as a function of the fiducial electronic density n°; solid line: 
for T = 3.71 x 10 6 K and Is, 2s, 2p, 3s, 3p, 3d states (from 
bottom to top); dashed line: for T = 2.32 x 10 5 K and the Is 
state. Upper and lower limits of energy bands are marked as 
e+ and e_. 



where w„ 



for the lowest level a of the ion with v 



electrons, and v* < v is the number of electrons remain- 
ing on the inner shells [y* = for Is and Is 2 ). 

The contribution of the internal degrees of freedom into 
the free energy is calculated as F; nt = — Y^ u N v ksT In Z v , 
where Z v = J2 a w va d va exp(— E va jk^T) is the internal 
partition function of the ion in the plasma. 



C. Total free energy model 

We evaluate F ee + Fa + F ie in Eq. (JIJ as described 
in Sec. Ill Al Albeit this is not strictly correct for ions 
with bound states, which are not pointlike, we need this 
approximation to make practical EOS calculations. 

The total free energy, normalized to Nik^T, can be 
written as 



— / + /int i 



Nik B T 

where f = £ x v f v + ft + Zf e 



fi - ln(raf) - 1 - 

fe=Xe~Pe, Xe 
-Pint 



Me 



n P knT 



/int 



(12a) 
(12b) 

(12c) 
, (12d) 

(12e) 



[to be used in Eq. Q] is Z v = Z n 



v*{\ 



and P e is the free-electron pressure. All terms of / can 
be calculated using the fitting formulae [25ll2r| . and only 
/int should be evaluated numerically. 
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III. THERMODYNAMIC EQUILIBRIUM 

A. Equilibrium conditions 

Thermodynamic equilibrium at constant V and T re- 
alizes at the minimum of the Hclmholtz free energy F tot . 
Since the total number of the ions in all ionization states 
is fixed, this minimum must be found under the con- 
straint N v — Ni. The charge neutrality condition 
is satisfied automatically, because the total number of 
electrons equals Z nuc in each ion cell by construction, 
however at cost of the a priori unconstrained volume. In 
order to maintain V = constant, one should impose the 
condition ^2^N u v u — V. 

These equilibrium conditions can be written as 

/ = minimum; Gn = Gy = 1; x v > 0, W; (13a) 
Gn = 2_, Xv ' Gv — ni^^x v v v , (13b) 

where / = f(xe,T, {x„}) is given by Eq. $T2$, m and T 
are fixed, while \e and x v may vary. 

B. Finding the equilibrium 

While solving the constrained minimization problem 
(|13fl . we take into account the condition Gjy = 1 explic- 
itly, by setting x — 1 — y^^" -1 x v , and discard those 
{x„} sets which would result in the negative right-hand 
side of this equation. As mentioned above, we do not con- 
sider the neutral atoms (v — Z nuc ). In order to satisfy 
the constraint Gy — 1, we use the Lagrange multiplier 
method. Namely, we minimize an auxiliary function 

$(n h T; X e, {xvh A) = / - AGy + \ 2 {G V - if (14) 

with respect to its arguments Xe and x v (1 < v < 
Znuc ~ 1) for different values of the Lagrange multiplier 
A, and find the A value that gives Gy — 1 at the mini- 
mum. The last (quadratic) term in Eq. i|14fl is an empir- 
ical regularization term which accelerates the solution. 
The solution provides the equilibrium values of .Ftot, Xe, 
x v , and n e = Zm. 

At each value of A, we approach min $ in two stages: 
first, a rough position of the minimum is found by the 
simplex method, and then it is refined by the Powell's 
conjugate-direction procedure [3^ |. In order to filter- 
out false local minima, the minimization procedure is re- 
peated several times with different initial sets of param- 
eters, and the absolute minimum is selected. A search 
for the root of the equation GV(A) = 1 is performed 
by bracketing and bisection, Because of the complicated 
dependence of / and Gy on the set of x„, and due to 
the limited accuracy of minimization, Gy(X) may ex- 
hibit a numerical discontinuity, which sometimes disal- 
lows the bisection to converge, so we have tried several 
initial guesses of A in such cases. 



C. Calculating thermodynamic functions 

Once ftot is calculated for a range of temperatures 
and densities, all thermodynamic functions can be found 
from its derivatives. The first derivatives give the pres- 
sure P = -dF tot /dV\ T , entropy S = -dF tot /dT\ v , and 
internal energy U = -T 2 d / dT {F tot / T)\ v = F+TS. The 
second derivatives give, for example, specific heat Gy = 
dU/dT\y and the pressure exponents (temperature and 
density logarithmic derivatives) xt = <9fnP/<91nT|y 
and Xp = —dhxP/d\n.V\T- In these derivatives, N{ is 
kept fixed, but x v depend on V and T, following the 
solution in Sec. lIIIBl 

Although such calculation looks simple, it is techni- 
cally complicated. We achieved the accuracy of / within 
0.003 over the p — T domain where the electron degener- 
acy is weak or moderate (x e S 10), and to four digits in 
the strongly degenerate regime (where Xe 3> 1), but this 
is insufficient for an accurate evaluation of the second and 
mixed derivatives of /. The difficulty is partly overcome 
by filtering the calculated values. We performed calcula- 
tions on a grid of (p, T) points and evaluated the deriva- 
tives at each (p, T) point from using the least-squares fit 
to the F values at a hundred of neighboring grid points. 

This filtering is not sufficient, if the electrons are 
strongly degenerate. In this case, the T-derivatives of 
In f tot, hit/, and InP are so small that a tiny numerical 
noise may preclude their evaluation. Fortunately, in this 
regime these derivatives are mainly determined by / fi . 
We use the following modification of Eq. I|12al) : 

/ = f + Smix + /', /' = J2 X » HXv/Zv)- (15) 

v 

The values of /', x v , and Z, and their p- and T- 
derivatives are calculated numerically, as described 
above, whereas / fi and its derivatives are obtained from 
the analytical fits |2fj. 

The calculated functions Cy, xt, and xp still exhibit 
a considerable numerical noise. To suppress it, we again 
employ the least-squares filtering. Improved values of 
pressure, consistent with the filtered xt, are obtained 
by numerical integration of the equation InP = J xt dt, 
starting from the lowest isotherm. 

The thermodynamic stability (Cy > 0, P > 0) and 
normality (xt > 0) require that S/Nik-g, monotonically 
increases with decreasing p or increasing T. To maintain 
these properties, we calculate S by integration of the 
equations dS/dV\ T = xtP/T and dS/dT\ v = C v /T, 
starting from the highest p and lowest T. 



IV. RESULTS FOR CARBON 

We have calculated the EOS for carbon at 2.34 x 
10 5 K < T < 3.63 x 10 6 K and 10~ 4 g cur 3 < p < 10 4 
g cm~ 3 ; at still higher p the carbon plasma is fully ion- 
ized. For lower temperatures, 7.5 x 10 4 K <T< 2.34 x 



6 




T (K) p (g cm" 3 ) 



FIG. 3: Number fractions of different carbon ions in the 
plasma (left vertical axis) versus temperature at p — 10~ 2 
g cm -3 . Solid lines: accurate results; dashed lines: Saha ap- 
proximation with current partition functions. Numbers of 
bound electrons v are marked near the curves. The thick 
solid curve shows the mean effective charge Z (right axis). 
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in Fig. [3] but at p = 1 g cm" 
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FIG. 5: Number fractions (left axis) and the mean effective 
charge Z (right axis) versus temperature at p = 100 g cm" 3 . 



FIG. 6: Same as in Fig. [£j but versus density at constant 
T = 10 6 K. 



10 5 K, we have calculated the EOS at 10 4 g cm 3 



< 



p < 0.04 g cm" 3 . 

Figures [2HS1 show the T-dependences of ion number 
fractions x v and the mean effective charge Z. At the 
lower densities, the electrons are nondegenerate. In this 
case, the mean ionization degree and Z depend sensi- 
tively on temperature. At the high density (Fig. 0, the 
electron degeneracy is significant, and the number of free 
electrons is mainly controlled by pressure, rather than 
temperature, so that Z varies weakly. However, the state 
of the bound electrons still depends appreciably on T: 
most of them are in the Is state at the higher T and in 
the Is 2 state at the lower T. 

The nonideality effects are less important at lower den- 
sity. Therefore, the abundance of individual ion species 
at low densities can be evaluated from the Saha equation 



Xu+l 



Z v 2 



m e K B 



3/2 



(16) 



as illustrated by Fig. [5] for p = 0.01 g cm -3 . Note that 
the shifts of bound-state levels in the plasma environment 
are included in Z v . At p -C 0.01 g cm" 3 , this approxima- 
tion and our calculations give identical results (this is one 
of the checks of our calculations), but at p > 0.01 g cm" 3 , 
Eq. (fT?)|l becomes progressively inaccurate (Fig. . The 
differences between Saha and our models in Figs. [3] and 
01 are due to the configurational effects (i.e., the devia- 
tions of the neutrality volumes from their ideal values; see 
Fig.nj and the Coulomb plasma nonideality (Sec. Ill A|) . 

The p dependences of the ionization states at T — 
10 6 K are shown in Fig. They exhibit pronounced 
maxima and minima due to the pressure ionization of 
particular bound states in particular ionization stages. 
These features are related to swelling and shrinking 
of the individual neutrality volumes relative to their 
rigid-background values (Fig. ^) and the corresponding 
changes in the internal partition functions. 

Figure \7\ presents normalized pressure as a function of 
density at T = 5 x 10 5 K. The vertical scale is smaller for 



FIG. 7: Isotherms of normalized pressure Pfriik^T for T — 
5 x 10 5 K. The present data (solid lines) are compared with the 
FGV77 (dot-dashed lines), FGV99 (dotted lines), and fully 
ionized plasma (dashed lines) models. Note the different scale 
in the figure above and below the horizontal long-dash line. 



FIG. 9: Temperature exponent \t for isotherms T = 2.4 x 
10 5 K and 1.26 x 10 6 K, compared with the FGV77 (dot- 
dashed lines), FGV99 (dotted lines), and fully ionized plasma 
(dashed lines) models (the curves marked with log T values) . 
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0.01 0.1 1 10 100 1000 
p (g cm" 3 ) 

FIG. 8: Isotherms of internal energy for T = 2.4 x 10 5 K, 
1.26 x 10 6 K, and 3.16 x 10 6 K (the curves are marked by 
log T values). The present data (solid lines) are compared 
with the FGV77 (dot-dashed lines), FGV99 (dotted lines), 
and fully ionized plasma (dashed lines) models. 



the upper part of the figure, to take account of the rapidly 
growing pressure of degenerate electrons. The difference 
between our results and the FGV99 tables is in general 
of the same magnitude as the difference between FGV77 
and FGV99. However, our isotherms exhibit more fea- 
tures. The slope of each isotherm varies near the den- 
sities where the ion composition of the plasma rapidly 
changes. These variations could not be revealed by the 
Thomas-Fermi model, but are easily grasped within the 
free-energy minimization technique. Related variations 
are seen in Fig. [SJ which shows isotherms of the internal 
energy per unit mass, U* = (U + Ua)/Nimi measured 
from the energy level —Uq of a nonionized ground-state 
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FIG. 10: Normalized specific heat at constant volume cv = 
Cv/Nik B for T = 3.16 x 10 s K, 10 6 K, and 3.16 x 10 6 K (the 
curves marked with log T values). Present data (solid lines) 
are compared with the fully ionized plasma model (dashed 
lines). 



carbon, which corresponds to a shift equal to 8.28 x 10 13 
erg g _1 with respect to the electron continuum level. The 
gap in the cold isotherm of the fully ionized plasma model 
(dashed line) corresponds to the region of instability of 
this model. Variations of the EOS due to the changing 
plasma composition with increasing density are also seen 
for the temperature derivative \t, shown in Fig.0 

Figure shows the normalized specific heat cy = 
Cy/NikB as a function of p for three values of T. At low 
p, where the electrons are nondegenerate, the ideal-gas 
value of cy is Cy — 1.5(Z + 1). Actually cy approaches 
Cy at T > 3 x 10 6 K, where the ionization is almost com- 
plete but in general, these two quantities can differ con- 
siderably, in particular because the heat absorbed by a 
partially ionized medium is spent not only on the increase 
of the kinetic energy of the ions and electrons, but also 
on thermal ionization. At very high densities, the specific 
heat is mostly determined by the ions. For a classical ion 
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solid cy = 3; for a strongly coupled ion liquid cy f» 3.4 
near the melting V (Eq. (17) of Ref. [2(|). The correc- 
tions due to the ion-electron interactions and quantizing 
ion motions are within 10% in the high-density part of 
Fig. This explains the limit cy ss 3 seen at high p. 

At p > 1000 g cm -3 , our model reproduces the EOS 
for a fully ionized nonideal plasma |2q . 

V. CONCLUSIONS 

We have developed a model for calculation of the EOS 
for dense, partially ionized plasmas, based on the free 
energy minimization method and suitable for pressure 
ionization zones. The free energy model is constructed 
in the framework of the chemical picture of plasmas and 
includes detailed self-consistent accounting of quantum 
states of partially ionized atoms in the plasma environ- 
ment. Occupation probability formalism based on the 
energy bands is used to treat pressure ionization. 

The developed technique is applied to the carbon 



plasma at temperatures T > 10 5 K, which is relevant 
for inner envelopes of the carbon-rich white dwarfs or 
outer accreted envelopes of the neutron stars. For gen- 
eral astrophysical applications of various type of stars, it 
is desirable to extend the calculated EOS to other chem- 
ical elements, first of all to oxygen. We are planning to 
perform such calculations in near future. 
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